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Abstract. The positions of atoms forming a carbon nanotube are usually described 
by using a system of generators of the symmetry group. Each atomic position 
corresponds to an element of the set Z x {0,1,..., n} x {0,1}, where n depends on 
the considered nanotube. We obtain an alternate rather different description by 
starting from a three-axes description of the honeycomb lattice. In our mathematical 
model, which is a factor space defined by an equivalence relation in the set 
Uvo, vi, v%) € Z 3 I vq + v\ + V2 & {0, 1}}, the neighbours of an atomic position can 
be described in a simpler way, and the mathematical objects with geometric or 
physical significance have a simpler and more symmetric form. We present some 
results concerning the linear representations of single-wall carbon nanotubes in order 
to illustrate the proposed approach. 
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1. Introduction 

A single-wall carbon nanotube is a cylindrical structure with a diameter of a few 
nanometers, periodic along its axis, which can be imagined as a rolled up honeycomb 
lattice. The high symmetry of carbon nanotubes has facilitated the theoretical 
investigation of the physical phenomena occurring in these materials P3 El EE EH El E] • 
The spatial symmetries (translations, rotations and screw axes, mirror and glide planes, 
etc.) form a line group, which is the maximal subgroup of the Euclidean group 
that leaves the nanotube invariant. The role of this group is analogous to that of 
crystallographic space groups in solid state physics. Some important properties of the 
band structure (electronic, phonon, etc.) can be directly deduced from the symmetry 
groups. 

The symmetry group of a carbon nanotube depends on the diameter of the tubule 
and on the helical arrangement of the carbon hexagons. The irreducible representations 
of these groups are well-known 001131, but, generally, equivalent representations offer 
distinct formal advantages. A calculation very simple in a representation can become 
much more complicated in an equivalent representation. Therefore, we think it is worth 
looking for new representations, and for new ways to describe the atomic structure of 
these materials. Our aim is to present an improved version of the mathematical model 
proposed in [3] and some applications illustrating this new approach. 

2. Honeycomb lattice in a three-axes description 

The vectors corresponding to the vertices of a regular triangle 

« = (!,<», -=(4f), -=(-i,-#) m 

allow us to define the bijecton 

£ — > L : (v , v 1 ,v 2 ) i-> v e + v\e\ + v 2 e 2 (2) 

from the set 

£={v = (vo, Vi, v 2 ) e Z 3 I v + v x + v 2 G {0; 1} } (3) 

to the set L of all the vertices of a honeycomb lattice. The subset £ of Z 3 becomes in 
this way a mathematical model for the honeycomb lattice. One can remark that 

£ = TU(T+ (1,0,0)) (4) 

where 

T = { v = (v , vi, v 2 ) e Z 3 | v + vi + v 2 = }. (5) 
The mapping 

d : £ x £ — > N, d(v, u) = \v — u \ + \vi — ui\ + \v 2 — u 2 \ (6) 
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is a distance on £, and u is a neighbour of order I of v if d(v , u) = I. The nearest 
neighbours of v are 

^° = (vo + e(v),vi,v 2 ) 

v 1 = {v Q ,v 1 +e(v),v 2 ) where e(v) = (-1)^+^+^ (7) 

v 2 = (v ,v 1 ,v 2 + e(v)) 

and the six next-to-nearest neighbours of v are the points i/- 7 = (V)- 7 corresponding to 
all the pairs with i ^ j. The symmetry group G of the honeycomb lattice coincides 
with the group of all the isometries of the metric space (£, d), and is generated by the 
transformations 

£ — ► £ : (v , v u v 2 ) i-» (v i, u 2 , u ) 

£ — > £ : (f , ui, v 2 ) h-> (u , t> 2 , ui) (8) 
£ — >£: (u ,ui,«2) •-»• {-Vo + -«2)- 
The subgroup of translations contained in G corresponds to T 

{ u | ve£ =>• t/+!iG£ } = { u G £ | e(«) = 1 } = T. (9) 

We can extend the description based on e , ei, e 2 to the whole plane. Each vector 
u admits the representation 

' 2 
v = ay^(v,ei)ei with a = - (10) 

O 

i=0 

and the usual scalar product and norm become 



) = a^2(v,ei}(u,ei}, 



{v,u i 

i=0 



2 

(I X 

i=0 

2 



The vectors eo, ei, e 2 form a tight frame and a system of coherent vectors j3] in 
For each vector v, the 'canonical coordinates' 

v = {v,e ), v 1 = (v,e 1 ), v 2 = (v,e 2 ) (12) 

satisfy the relation v + v\ + v 2 = 0, and we can identify M 2 with the space 

K = {k = (k , h, k 2 ) | fco, h, k 2 eR, k + k x + k 2 = } (13) 



by using the linear isomorphism 



2 



/C — >R 2 : (ko,h,k 2 ) ^aJ2 k i e i- ( u ) 

i=0 

The representation of a vector v as a linear combination of eo, ei, e 2 is not unique 

2 2 

v = a -Oj €i = a ^^(wj + a) e i (15) 

i=0 i=0 

for any a E R. All the elements (v , t>i, v 2 ) of the set 
{(v + a, -Oi + a, v 2 + a) a G K } 
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correspond to the same point of the plane, and 

2 2 2 

(k, v ) = a kiVi = a k^Vi + a) = a k(Vi (16) 

i=0 i=0 i=0 

for all k G /C. 

3. An alternate mathematical model for carbon nanotubes 

A single-wall carbon nanotube can be visualized as the structure obtained by rolling a 
honeycomb lattice (which is a mathematical model for a graphene sheet) such that the 
endpoints of a translation vector c G T are folded one onto the other (see figure 1). 
The vector c is called the chirality of the tubule. After the graphene sheet rolling, the 
points v — 2c, v — c, v, v + c, v + 2c, ... are folded one onto the other, for any v G £. 
Therefore, each element of the set [v] — v + Zc, that is, each element of the set 

v x , v 2 \ = { (v , v 1: v 2 ) + j(c , ci, c 2 ) I j G Z } (17) 

describe the same point of the nanotube. The subset of the factor space Z 3 /Zc 

Z 3 



, + , 1 + ,2G { 0;l}| (18) 

can be regarded as a mathematical model for the nanotube of chirality c. A symmetry 
transformation C — > C : v h- > (?f of the honeycomb lattice satisfying the relation 

[ v ] = [«] =► M = H (19) 

defines a symmetry transformation of nanotube £ c , namely 

j:4-»£ c :[i;]h [^]. (20) 

Each nanotube C c admits the symmetry transformations: 

t : C c — ► C c : [w , ui, u 2 ] >-> [-uo + 1, -Vi, -v 2 ] ^ 
g u : C c — >■ C c : [v] h- > [v + u] for any u E T. 

Let n be the greatest common divisor of Co, ci, c 2 , and let c = (co, ci, £2), where 

111 

c = -c , ci = -ci, c 2 = -c 2 . (22) 
n n n 

The transformation 515 represents a rotation of nanotube of angle 27r/n with respect to 

its axis. Since 

(c 2 - ci)c + (c - c 2 )ci + (cx - c )c 2 = (23) 

the vector (c 2 — ci, Co — c 2 , ci — Co) is orthogonal to c (it has the direction of nanotube 
axis). The transformation g t corresponding to 

* = i(c2 - cx, c - c 2 , cx - c ) (24) 

where 

jf g 2 Cx G 3Z 

ft = gcd{c 2 - ci, c - c 2 , ci - c } = <( . 2 } (25) 

it c 2 — cx f. 3Z 
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Figure 1. The unit cell of the carbon nanotube with the chiral vector c = (10, —2, —8). 
In this case n — 2, c = (5,-1, -4), t = (— 1, 3, -2), to = (1,0, -1) and q = 14. 



represents the shortest pure translation of nanotube. From c + C\ + c 2 = we get 

(ci - c 2 ) 2 + ( C2 - Co ) 2 + ( Co - Cl ) 2 = 3(c 2 + cj + cj) (26) 

that is, 72 2 ||t|| 2 = 3||c|| 2 , whence 

q=^-(cl + cl + cl)enZ. (27) 

For any u G T the projections of u on c and £ can be written 

(u,c) ( Cx-Cq c 2 -c \c (u,t) ( c 2 d\t 

" c = "i— ^ — + M 2^5 — -, -f-— -t = f «i w 2 — J- (28) 



||c|| 2 \ 72. 71 J q \\t\\ 2 \ n nJ q 

where q — -. Since 

gcd{(ci - Cq)/H, (c 2 - c )/72} = 1 and gcd{c 2 /n, Cx/n} = 1 (29) 

it follows that the projection of T on c is Z-, and the projection of T on t is Z4. Let 
wGTbe the shortest vector with 

(W ' t} t=-. (30) 



From the relation 

» = feF c+ W' (31) 

we get [gw] = [nt]. 

Without lost of generality we can assume c > C\ > c 2 . In the case Ci = c 2 we have 
an armchair nanotube, and in the case Cx = a zig-zag nanotube. The nanotubes with 
7^ ci 7^ c 2 are called chiral nanotubes. Our approach works for any single-wall carbon 
nanotube, but in this paper we restrict us to chiral nanotubes. The symmetry group 
G c of the nanotube C c is generated by the transformations g — o — g w and r 

G c = (g, a,r \ go = ag, g n = r 2 = (ar) 2 = (gr) 2 = e). (32) 
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For each [v] G C c there exist s G Z, m G {0, 1, n — 1} and p G {0; 1} uniquely 
determined such that [v] = a s g m r p [0, 0, 0]. The usual description 0113 of the atomic 
positions of the atoms forming a carbon nanotube is based on this remark, and the set 

{(s,m,p) | s G Z, m G {0, 1, n - 1}, p G {0,1}} (33) 

is used as a mathematical model. The alternate mathematical model we present in this 
paper is £ c . 

4. A tight-binding approach to carbon nanotubes 

Consider the Hilbert space (7 2 (£ c ), (, )), where 



Ei'^)i 2<o ° ( 34 ) 



v€C c 



<^2> = X>i00^2(tO (35) 

and the unitary representation of G c in / 2 (£ c ) defined by 

5 : l 2 (C c ) — > Z 2 (£ c ) (gVOM = V'Gr'M)- (36) 
For each « G (0, oo), the linear operator 

2 

if : l 2 (C c ) — / 2 (£ c ) [«] = « J] (37) 

3=0 

is a self-adjoint G c -invariant operator. 

If k G /C is such that (fc, c) G 27rZ then 

e i<M = e Hk,v+jc) for any ^ G z (gg) 

and hence, the function 

C c — > C : [«].-> e -1 ^^ (39) 

is well-defined (it does not depend on the representative we choose for [v]). The 
Hamiltonian used in the tight-binding description of tc bands has the form ([37)1 . 

Theorem 1. a) For any k such that (fc,c)G27rZ the numbers ±E(k), where 

E(k) = K,\e ikoa + e ikia + e ik2CL \ 

. (40) 

= Ky3 + 2 cos(A;o — k\)a + 2 cos(A;i — k 2 )a + 2 cos(/c2 — ko)a 

belong to the spectrum of H . 

b) The bounded functions belonging to an extension ofl 2 (C c ) 

— C, = e-^VfM (41) 

where 

V±\v]-{ lf £{V) = 1 (42) 
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( - i are ( e ifc ° a + e ifcia + e ik2a ) % f e ik,a + e ikia + e ik2a ^ 
X(k) = \ 2 V ' (43) 

[ if e ik ° a + e ifcia + e ik2a = 

are eig en functions of H corresponding to the eigenvalues ±E(k), that is, 

Hij,£ = ±E(k)ift. (44) 

Proof, a) The function ip : C c — > C, ip[v] = e -1 ^' v ^[i;], where 

r -1 a if e(v) = 1 , 
*> = { /> if e(J) = -1 (45) 

and a, f3 are two constants, satisfies the relation Hip = Eip if and only if (a,/3) is a 
solution of the system of equations 

K(e~ ikoa + e~ ikia + e- ik * a ){3 = Ea 

K(e ikoa + e ikia + e ik2a )a = E(3. ^ ' 



This system has non-trivial solutions if and only if 

-\k\a _|_ p — \k2a\ 

(47) 



—E K(e~ ikoa + e~ ikia + e~ ik2a ) 

n(e ik ° a + e ikia + e ik2a ) -E 



that is, if and only if E is one of the numbers ±E(k). 

b) If e ife ° a + e ifcia + e ifc ° a ^ then the equation K(e ik ° a + e ikia + e ik2a )a = ±E(k)f3 leads to 



gifcoa _|_ gifcia _|_ gifoa 



P = ±a—r — '—r — '—r— = ±a e" 2iA(fc) . (48) 

gifeoa _|_ gifeia _|_ gife2a ' 

Choosing a = e lX ^ we get (3 = ±e~ lX ( k \ and hence, up to a multiplicative constant, the 
solution of Hi/j = ±E(k) ip is ip = i/j£. If e ifc ° a + e ifcia + e ifc ° a = then E(k) = 0, and the 
functions 

are eigenfunctions. ■ 

The relation (k,c) G 27rZ defines a family of equidistant straight lines orthogonal to c 
with the distance between neighbouring lines equal to 27r/||c||. 

Theorem 2. a) The function 

E-.K — ►[0,3k] E(k) = K \e ikoa + e ikia + e ik2a \ (50) 
is even and periodic 

E(k) = E(-k) E(k) = E(k + b ) = E(k + b l ) = E(k + b 2 ) (51) 

where 

/4vr 2tt 2tt\ / 2tt 4tt 2tt\ / 2vr 2tt 4tt\ 

\3a' 3a' 3a/' ' \ 3a '3a' 3a/' 2 \ 3a' 3a' 3a/ 
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Figure 2. The hexagonal domain B, the points of A (indicated by •) and the points 
k with E(k) = 0, called K points (indicated by o). 



b) The functions ip k — i^ k and ip k are eigenfunctions of g and a 

g = J<M a = (53) 

satisfy the relation 

Ti>± = ±e- ik ° a i>t k . (54) 

c) The eigenspaces 

£ k = { a^ k + (3^)- k | a, (3 e C } , . 

£- = {ar k +^Z k I a, PeC } [ } 

corresponding to E{k) and —E(k) are G c -invariant, £ k = £- k and £ k = £Z k - 

d) If k is such that E{k) ^ then 

Vf +6j = e i7r/3 ^ £ k+bi = £ k (56) 

for any i G {0, 1, 2}. 



Proof. Since e l47r / 3 = e l27r / 3 we have 

E(k + b ) = K \e ikoa e i4n/3 + e ikia e~ i2n/3 + e ik2a e~ i2n/3 \ = E{k). 

X(k + bo) = -\ arg [(e ife ° a + e ifcia + e ifc2a ) e" i2w/3 ] = X(k) + \ 

for any k with E(k) ^ 0. Iff is such that e(v) = 1 then (b ,v) = 2ttv and 

il> k+bo [v] = e - i{k+bo ' v) e lX{k+bo) = e- i{k ' v) e" i2 ™° e iA(fc) e m/z = e m/z ip k [v}. 

If v is such that e(v) = —1 then (b , v) = 2-kvq — 2ir/3 and 

il> k+bo [v] = e - i{k+bo ' v) e ~ iA(fe+6o) = e~ i{k ' v) e" 127 ™ e i27r / 3 e - iA(fc) e" i7r/3 = e m/z ^ k [v\. 
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From the periodicity of E(k) and £ k it follows that we can restrict our analysis to the 
case k G B c: where 

B c = { k G B | (k, c) G 2nZ } (57) 

and B is the hexagonal domain (see figure 2) 



B=lke)C 



2ir 2-n 1 

■ ^ < h < — for any i G {0, 1,2} L (58) 

Theorem 3. For even n the space £ k is one- dimensional if and only if k G A, where 

a — I (o o o) ± ( 27r 71 n \ ± ( 71 27r ± ( 71 71 ^\ 

\ ' ' ' \ 3a' 3a' 3a/' \3a' 3a' 3a/' \3a' 3a' 3a/ 
Proof. If £ k is one-dimensional then there is a constant C such that 

= Cipk[v] for any v G £ c . 

For v = (0, 0, 0) we get C = e~ 2iX( - k \ Therefore 



c 



2i<fc,«> = 1 if £ t v \ = x 



^ = C ^ ^ < e H^ =C 2 if £ (v) = -l. (59) 

Particularly, we must have e 2lk ° a = e 2lkia = e 2lk2Cl . For k e B this relation is possible 
only if there are a,/? 6 {0, ±1} such that k\ = k + an /a and /c 2 = k + fin/a. From 
+ k\ + A; 2 = we get 

/ a + p7i 2a- pn -a + 2pn\ 

k = 5 , — 5 , 5 • 60 

\ 3a 3a 3 a/ 

All these points points lie on the family of straight line defined by (k, c) G 27rZ since 
(k, c) = koCoa + (k + a—^j c\a + (k§ + c 2 a = (aci + {3c<i)tx G 2ix7h. (61) 
In the case a = {3 = we get k = (0, 0, 0), ip k [ v ] — 1> an d 

e^fc = i>k o-ipk = i) k ri) k = ip k . (62) 

In the case a = (3 = 1 we get k = (— 27r/3a, 7r/3a, 7r/3a), V'fcM — ( — l) v °e~ 17r ^ 6 , and 

£^ fc = (-l) £( Vfc ^ = (-1)"°^ ^ = -i, k . (63) 

All the other points k G B c with dim^ = 1 can be obtained from the analyzed points 
by using permutations of coordinates and/or a multiplication by (-1). ■ 

In the case when n is odd only a part of the points of A lies on the family of straight 
lines defined by (k, c) G 2nZ,. 
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5. Two-dimensional representations of G c 

If k G B C \A then dim^ = 2 and the matrices of g, a and r in the basis {tpk, ip-k] are 

\ / e i{k,w) q \ / q e ik a \ 

e -i(k,c) y ° - I q e -i(fc,»> 1 ' T - l e -ifcoa j ^ ' 

Theorem 4. For eac/i k G £> C \A t/ie representations of G c in £ k and £ k are equivalent. 
They are reducible if and only if the numbers (k,c) and (k,w) belong to Jm. 

Proof. The linear transformation — ► £ k : aipk + l— * i^V^ ~~ ift^-k * s an 
isomorphism of representations. If 




P 



a P 
7 5 




e -i(k,c) J I g e -i<fc,c) 



e i{k,w) g \ / e i(fc,u>) g 






is the matrix in the basis {ip k , V'-fc} °f the projector corresponding to a G c -invariant 
subspace of then P 2 = P and 

n \ / „i(k,c) n \ / „ « \ 

(65) 
(66) 
(67) 

\ e " u I \ 1 I 
From the last three relations it follows that 

pe i{k ' £) = pe~ i{k ' £) pe i{k ' w) = pe~ [{k ' w) ae [k ° a = 5e ik ° a (68) 
^ e i(fc,c) _ ^ e -i(fc,c} ^ e i{k,w) _ ^ e -i(k,w) [3e~ ikoa = 7e ifc ° a . (69) 

Since e lfc ° a ^ we obtain a = 6. If either (k,c) G" Z7r or (k,w) G' Z7r then (3 = 7 = 
and the representation is irreducible. If the numbers (fc, c) and (fc, if) belong to Zir 
then imposing the condition P 2 = P we obtain 

/ x e ifc a \ 1/1 -e ifcoa \ 

P = - orP = - . (70) 

2 \ g-ifcoa J / 2 \ _g-ifcoa ]_ / 

These complementary projectors correspond to the decomposition 

= MVk + e- ifco >_ fc ) I a G C } © {a# fc - e- ifc °>_ fc ) | a G C } (71) 
of £k into direct sum of one- dimensional G c -invariant subspaces, and 

g(ij k ± e- ifc °>_ fc ) = (-l) m (^ fc ± e- ifc °^_ fc ) 

<r(V>* ± e- ifco >_ fe ) = (-lWfc ± e- ifcoa ^-fc) 

r(^ ± e- ifc °^-fc) = ±(^A ± e" ifc0 ^_ fc ) 
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where m, p G Z are such that (k, c) = ran and (k, w) = pn . ■ 
The relations (JMj) define for each k belonging to the set 

B vrred = { fc g | ^ ^ £ ^ Qr w ) £ Z 7T } (72) 

a two-dimensional irreducible representation T> c (k). Some of these representations are 
equivalent. Particularly, V c {k) = V c (—k). 



6. Clebsch-Gordan coefficients 

Let k, k! G i% rred be such that k + = k + k' and k~ = k - k! belong to £f rec! . The direct 
product of the representations T> c (k) and T> c {k!) admits the decomposition 

V c (k) ® D c (fc') = X? C (A; + ) © D C (AT) 

and the matrices corresponding to q, o~ and r are 

\ 



j(k,c) Q 



3 i(fc',c) Q 







-i(fc,c> 







-i(fe',c) 







-i(A;,to) 







-i{fc',w> 



e 

e -ifcoa Q 



e 

e -ifc'oa g 



e i(fc+,c> 















e 1 


(fc- 















e~ 


i(fc",c) 















/ s k+ > 


w) 















c 1 


(k~,w) 















e -i(fc- 


V o 























pi/CQ"(2 















e" 




a 








ik^a 



















o \ 






(73) 



respectively. The unitary matrix 

/ 1 \ 
10 



M 



1 
\ 1 / 



(74) 



satisfies the relations 

/ e i < fc+ ' a > 



M 



-i 



v 


V o 



i(k ,c) Q 

o e -'^ k ~^ 





o \ 





g-i(fc+,c) J 



/ 



M 






V o 



\{k+,c) Q 



o \ 









3 i(fc ,c> 







i(fc-,c> y 
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o \ 








(k+,w) 












-l 







e i(k~,w) 





















g-i(fc+,uj) 







M~ 
















M= 


























e -i<*- 


,w) 



















e i(fc- 


,w) 






\ 













e" 1 


{k+,w) j 






V 

















( 













a \ 




( 













\ 




-l 



















e" 


-ik^a 













M~ 














M = 
























































V 


e -ifeJo 











y 1 




V 













J 
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\ 



-i(fc ,w) i 



(75) 



Therefore, the entries of M are Clebsch- Gordon coefficients jTTj corresponding to the 
considered direct product. In this case, the only non-null coefficients are 

{kk'll\k + l) = (kk'12\k-l) = (kk'21\k-2) = (kk'22\k + 2) = 1. (76) 

More details concerning the Clebsch-Gordon coefficients and their applications in carbon 
nanotube physics can be found in the articles of Damnjanovic et al. 17| 

7. Concluding remarks 

The present paper can be regarded as a pure mathematical exercise. We have defined 
the factor sets £ c , the groups G c acting on C c as groups of permutations, and we have 
studied certain representations of these groups defined on some spaces of functions 
if) : C c — > C. We have proved that the groups G c are isomorphic to the symmetry 
groups of single-wall carbon nanotubes, and the considered representations are directly 
related to some representations used in carbon nanotube physics. 

The present paper can also be regarded as presenting an alternate mathematical 
model for carbon nanotubes. We think that this alternate approach offers some formal 
advantages, namely, certain calculations may be simpler in this approach than in the 
usual one. 
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